Efimov States of Heavy Impurities in a Bose-Einstein Condensate 



(N 

O 

<N 

o 



C/3 

a- 



O 

> 

O 
(N 



X: 



N. T. Zinner 

Department of Physics and Astronomy, Aarhus University, Aarhus C, DK-8000 
(Dated: December 10, 2012) 

We consider the problem of two heavy impurity particles embedded in a gas of weakly-interacting 
light mass bosonic particles in the condensed state. Using the Bogoliubov approach to describe the 
bosonic gas and the Born-Oppenheimer approximation for the three-body dynamics, we calculate 
the modification to the heavy-heavy two-body potential due to the presence of the condensate. For 
the case of resonant interaction between the light bosons and the impurities, we present (semi)- 
analytical results for the potential in the limit of a large condensate coherence length. In particular, 
we find a formula for the modification to the Efimov scaling factor due to the presence of a degenerate 
bosonic gas background. 

PACS numbers: 03.65.Ge,03.75.Hh,67.85.Pq 



I. INTRODUCTION 

The three-body spectrum for two heavy and one light 
particle with short-range resonant interactions between 
the heavy and light particles can be treated using the 
Born-Oppenheimer approximation. One decouples the 
slow degrees of freedom (relative motion of the two heavy 
particles) from the motion of the light particle to obtain 
an effective heavy-heavy potential. What is found is the 
famous Efimov spectrum of geometrically scaling three- 
body states where subsequent states differ by a simple 
factor of the form e^/^° where sq is called the scaling pa- 
rameter For equal mass particles, e^/^° ^ 22.7, which 
is a severe limiting factor in the quest for observing many 
states. However, in the heavy-heavy-light system the fac- 
tor can be much smaller 0, Q ■ This has generated a lot 
of interest in studying three-body systems with mixed 
atomic species in the ultracold regime 

The ultracold gases used for three-body experiments 
are typically at very low densities such that the three- 
body dynamics can be described by assuming that three- 
body states are essentially isolated entities in spite of the 
fact that there is always some finite density background 
present. However, experiments should be able to push 
into higher densities and thus be able to investigate the 
effects of the finite background on few-body properties. 
Some recent works have been devoted to the question of 
what happens when the particles are fermions that have 
to obey the Pauli principle . Here we will address the 
related issue of what happens when there is a Bose gas 
background that is in the condensate phase. A few other 
papers have recently discussed this type of question from 
the many-body point-of-view Here we take a 

different point-of-view and consider the background con- 
densate as a perturbation on the three-body state. For 
simplicity we will consider the case of two heavy impurity 
particles in a Bose gas of light bosons. For the two-body 
case the problem of an impurity in a condensate with 
strong interactions was revived recently fl34T^ , and like- 
wise for the case of multiple such impurities (I5i, ,16|. The 
Bose-Fermi mixture in cold atomic gases was discussed 
in detail about a decade ago [13, [3 and the later works 



address a particularly interesting regime with large im- 
balance in the populations. However, we are not aware 
of any works discussing the presence of an Efimov effect 
in a condensate background and the influence that the 
background can have on the universal three-body spec- 
trum. 

Our approach will be based on the Bogoliubov approx- 
imation for the condensate dynamics. We will work in 
the limit of a weakly-coupled Bose gas with repulsive 
short-range interactions, i.e. small and positive scatter- 
ing length, aB > 0. The impurity-light boson interaction 
will be allowed to take on any value. The condensate 
coherence length, ^ — 1/ ^/ETm^OB with riQ the conden- 
sate gas density, will therefore be large and its inverse 
a useful expansion parameter. The problem we solve is 
very similar to the classical problem of electron-electron 
interaction mediated by phonons; the impurities are the 
(heavy) electrons while the light bosons are the phonons. 
Of course, in the case of a Bose condensate of massive 
particle, the phonon dispersion is only linear for low 
momenta and eventually becomes quadratic. This will 
be properly taken into account in our framework. The 
problem considered here is closely related to the bosonic 
Kondo problem with two impurities and a bosonic bath 
P^I - |23 |. The study here should thus be of interest for 
both cold atomic gases and for condensed-matter physics. 



II. THEORETICAL MODEL 

We consider two heavy impurities of mass M and a gas 
of light bosonic particles of mass m ^ M . The Hamilto- 
nian is 

+Ub Eq nB{q)nB{-q) + Uib Eq nB{q)ni{-q), (1) 

where c are the impurity and h the boson operators. 
We use zero-range density-density impurity-boson and 
boson-boson interaction with nj{q) = Efe '^I+q'^fc ^^d 
"5(9) = Efc^I+o^fc- The dispersions are e/(fc) — 



frk^ /2M and es(fc) — hrk^ /2m. In the weakly-coupled 



2 



limit we have Ub = 47r S^as/m. The paramter Uib will 
be discussed later. We use Bogoliubov theory to describe 
the light bosonic particles igSj. We therefore transform 
to quasi-particles, 7k and 7^, in the standard way. Fur- 
thermore, we will assume that the condensate density is 
small so that the number of quasi-particles is also small. 
This allows us to neglect all terms expect 7^7fc in the 
transformed Hamiltonian. Dropping unimportant con- 
stant terms, the bosonic dispersion and the interaction 
term becomes 
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where Esik) = ^UBnoh?k^ /niB + {h?k^/2mBY. This 
corresponds to impurity particles with dispersion ei{k) 
interacting with Bose gas particles with dispersion E{k) 
through a contact interaction with strength Uib- In 
the case where E{k) is linear in fc, this corresponds to 
a system of (heavy) electrons interacting with phonons 
through a non-dispersive zero-range intearction. We now 
proceed to solve the three-body problem in this setting. 



A. Born-Oppenheimer Approximation 

As we have just argued, we can model the in- 
teraction of the impurities and the bosons via a 
zero-range interaction, and we thus write V{r) = 
Vo [6{r — i2/2) + 5{r + ii/2)], where we assume that the 
two heavy impurities are located at ±R/2. This poten- 
tial needs to be regularized since as it stands it leads to 
an ultraviolet divergence. We return to this point below. 
The essence of the Born-Oppenheimer approximation is 
that we first solve for the dynamics of the light bosonic 
particles while assuming that R is fixed, and then pro- 
ceed to consider the Schrodinger equation for the two 
heavy particles as function of R. The relative distance 
of the heavy particles, R, is our adiabatic variable which 
we assume changes on a much slower timescale than the 
positions of the bosonic particles. 

Now we consider the Schrodinger equation for the par- 
ticle of mass m in this potential 



(3) 



where E^ is the energy and the wave function of a 
light bosonic (quasi)-particle. Since V{r) contains delta- 
functions, it is convenient to work in momentum-space. 
The Schrodinger equation in momentum-space can then 
be written 

E{k)cf,{k) + j^ J d^k'<j>{k')V{k-k')=ERcl>{k), (4) 

where V{q) — 2Vq cos (g ■ y) • Integration over fc on both 
sides, assuming that (j) is an even function of k (s-wave 
solutions), and elementary trigonometric manipulations 



reduces the equation to the form 
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E{k)-ER (27r)' 



cos (fc • R) 



E{k) - Er ■ 
(5) 

The time has now come to relate Vb and the scattering 
length of the interaction between the impurity and the 
bosonic particles, a. While this can be done most ele- 
gantly by using Tan's pseudopotential [1^, 113], we use 
a traditional approach that is very easy in the Born- 
Oppenheimer limit. From the Lippmann-Schwinger 
equation for the impurity-boson scattering we have 



1 



1 
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(6) 



where ^ = mM/{m + M) and = h^k^ /2fi are reduced 
mass and energy. Since we assume m <^ M , we can safely 



use = Efe and 11 



Note that we do use the bare 



dispersion of the bosons, h^k^ /2m, and not E{k). This 
is necessary since the heavy-light scattering length, a, is 
defined in vacuum and the Lippmann-Schwinger prob- 
lem must therefore also be solved in vacuum. Inserting 
the relation between Vo and a we arrive at our central 
equation 
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[x4+A2x2] 1/2 + 1' 
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where we have defined = —2mER/h^ and A — l/{aS,). 

In the case where the bosons are non-interacting, i.e. 
as — ?> 0, we have ^ — >■ c» and thus A — >■ 0. In this limit 
the integrals can be performed analytically and we arrive 
at the well-known formula 
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aR = - + e-"-^. 



(8) 



In the case of resonant interactions, \a\ = 00, the equa- 
tion has the solution aR — xq ^ 0.567. To lowest order 
in R/a (i.e. for large \a\ and/or small radii), we find 



2^2 r 



Ef 



K'xl 
2mi?2 



1 



1 



R 



xf^{l + e^o) a 



(9) 



Another interesting limit, is i? 3> a. In that case we can 
neglect the exponential term in Eq. and we find a 
well-known result 



Er, 



2ma2 



(10) 



Notice that this only works for a > 0, since the a < 
case has no solution. This energy is the usual energy of 
a particle of mass m in the delta-function potential of 
a much heavier particle of mass M ^ m (i.e. a fixed 
potential center). The physical interpretation is that the 
small mass particle forms a bound state with one of the 
heavy particles. 
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FIG. 1: Solution of the heavy-light two-body problem in the 
presence of a condensate of light particles. The solid (blue) 
line shows the energy in the presence of a condensate, while 
the dashed (black) line is the corresponding vacuum solution. 
For I S> 1, the condensate and vacuum solutions become 
indistinguishable. The dotted (red) curve is obtained from 
Eq. (|12p . The atom-dimer continuum is on the right side 
of the curves, while we expect universal three-body states 
to appear on the left. A clear change of threshold at small 
binding energy is found due to the condensate background. 



III. HEAVY-LIGHT TWO-BODY STATES 

We first consider a single impurity. This problem can 
be handled in the same manner as the discussion above. 
The only difference is that in Eq. ([7]) the integral with 
the sine term is absent. The first integral in Eq. (O is 
analytically tractable but the expression is long and cum- 
bersome and does not really yield any insights. However, 
since we are concerned with the large ^ limit we can make 
an expansion in A inside the integral. This requires the 
stronger condition of large, which must be checked 
after doing the integral. After expansion one has 



dx 



- 1 



(11) 



and using this we arrive at i = a+-r^ which is a hidden 
second degree equation in a with one positive root 
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(12) 



which requires ^/a > 1. 

In Fig. [1] we show the two-body binding energies of the 
heavy-light system with and without the presence of a 
condensate. We see a clear tendency of the condensate 
to push the threshold away from unitarity and into the 
regime of small positive a, i.e. one needs stronger attrac- 
tion to bind the light-heavy system in the presence of 



a condensate. A shift of the two-body threshold is also 
found with fermionic backgrounds However, in the 

light fermion and heavy impurity case the effect is op- 
posite when the Fermi sea is inert (with no particle-hole 
excitations) [Tj] and pushes the threshold to the a < 
region, i.e it enhances binding. This can be understood 
in the following way. In the case of a Fermi gas and an 
impurity, the binding is provided by the fact that there is 
a Fermi surface (reducing a three-dimensional problem to 
a two-dimensional one where it is much easier to produce 
bound states). Around the Fermi surface we have a linear 
dispersion of particles. This is the same as the situation 
here with bosons that have a linear dispersion around 
zero momentum. The Fermi surface moves closer to zero 
momentum as the mass of the fermions increase. So we 
see that there is a mapping to the problem with a Fermi 
sea instead of a condensate, but light bosons map to heavy 
fermions. Whether this mapping holds more generally 
away from the extreme light/heavy mass regions cannot 
be addressed within the current Born-Oppenheimer ap- 
proximation. This is an interesting question for future 
studies. 



IV. TWO IMPURITIES IN A CONDENSATE 

We now proceed to consider the effect of a condensate 
of light particles on two impurities. The solution for the 
potential of the two heavy impurities, given implicitly 
through a, is found by solving Eq. ([7]). While the first 
integral in Eq. ([7]) is exactly solvable, the second term is 
very difficult to handle both analytically and numerically. 
We have solved it using various approximation schemes 
and present the results in Fig. [2l The lines termed 'full 
solution' are expansion of Eq. ([T]) to 2nd and 4th order in 
A. The full solution to second order yields the equation 



aR + -aRA^ = exp(-ai? - -aRA^) 
4 4 ' 



(13) 



from which the 'lowest Taylor' is obtained by second or- 
der expansion around i?/^ = 0. In the absence of the con- 
densate, the solution is given by the first term in Eq. ^ 
which can be written aR = xq. We thus see that the 
effect of the condensate and the corresponding change 
in the dispersion of the light particle, is to suppress the 
attraction between the impurities. The scale of this sup- 
pression is not surprisingly ^. We arrive at the conclusion 
that the presence of a condensate of light particles will 
tend to make universal three-body states more difficult 
to form in this non-trivial many-body background. 

Numerically, it becomes increasingly difficult to find 
a solution as i?/^ — >■ xo, and beyond this point we have 
failed to find solutions. This indicates that the expansion 
of the integrand before performing the integrals becomes 
problematic. Fig. [5] does indicates that the quantity aR 
goes to zero at some finite value of i?/^. Numerically, we 
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FIG. 2: aR as a function R/S,- A full solution could only be 
found for i?/^ < xq- The full (blue) line is the full solution to 
2nd order in A, while the dashed (red) shows the full solution 
to 4th order. For comparison, the dotted (black) line shows 
the second order expansion around zero. The thin dotted lines 
guide the eye to the positions beyond which no full solutions 
were found. 



find that a good fit to the solutions for aR > xq /2 is 

1/2 



aR = Xq 



1 



R 



(14) 



Another way to approach the limit of where becomes 
small is to consider the case ^ 1 and make this ap- 
proximation before doing the integrals in Eq. [T] The first 
integral is analytic and can be easily expanded, while the 
second integral can be handled by using ^/x'^ + A^x"^ — 
Ax for 3> 1, taking a sine transform and subsequently 
doing the expansion. This yields the expression 



1 
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(15) 



which is always negative. This indicates that the Er 
does indeed change sign as i? ^ ^. However, taking the 
limit of i? — ^ in this result yields — ^ and A^ is now of 
order one which is in conflict with the limit we used to 
derive the expression. However, it strongly suggests that 
En will go to zero for some value of i?/^. 

Our results for the potential between the two heavy 
impurities, Eji, are plotted in Fig. [3l We see that the 
presence of the condensate tends to make the potential 
go to zero faster than R^^ as R increases. An effective 
repulsive effect is thus seen to originate from the many- 
body background in the light particle component. The 
Efimov effect depends on the R~^ functional form of the 
potential, and of course on the fact that it is attractive 
[3, Q . Our results imply that the condensate coherence 



FIG. 3: Effective potential for the two heavy impurities, Er, 
as a function of their relative distance R. The solid (blue) 
curve shows the result from a second order expansion of 
Eq. with respect to (a^)~^, while the dashed (black) curve 

2 

gives the results with no condensate background ( — ■^). The 
end point of the solid curve is at R/^ = xo where the second 
order solution ceases to exist. However, our results generally 
indicate that the potential El goes to zero faster than R~^ 
for large R/C 



length, ^, must be considered when estimating the Efi- 
mov effect in a degenerate Bose gas setting. Moreover, 
the length scale of the modification is ^ to within a fac- 
tor or order one. We therefore conclude that at unitarity, 
the number of universal three-body bound state can be 
estimate by analogy with Efimov's original formula [l[ 
and becomes 



N, « f^^log ( ^ 
Ho 



(16) 



where sn is the scale factor and Rq is a short-distance 
cut-off [1]. More generally, we expect that whichever is 
smaller of \a\ and ^ will cut off the number of states 
allowed in the spectrum. 

We note that our findings do not indicate any modifica- 
tion of the potential at small distances. This can be un- 
derstood from the dispersion relation of the light particles 
which become quadratic at high momenta corresponding 
to short distance. Condensates modify behavior at long- 
range (where low-energy universal three-body states may 
reside), short-distance is left more or less undisturbed. 
This means that the short distance cut-off that is neces- 
sary to bound the attractive R~^ and avoid collapse is 
provided by short-distance physics (i?o in Eq. (|16|) . of- 
ten called the three-body parameter). However, recent 
studies have found a strong connection between Rq and 
the two-body physics of the systems as given by the van 
der Waals length, rvdw |28l - l34| . In the case studied here 
where the impurity is much heavier than the light con- 
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stituent, Rq is given simply by the details of the heavy- 
heavy system as is seen from the derivation of Eji above 
(and from the numerical results of Ref. [35|). We there- 
fore expect that Rq ~ fvdw to within factors of order 
unity (between [s^ and 2 [s^). Note here that we 
are assuming that there are no resonances in the heavy- 
heavy system. This should always be closely investigated 
in realistic setups where interactions are typically con- 
trolled by external fields [s^ that affect both heavy-light 
and heavy-heavy systems. 



V. EXPERIMENTAL CONSIDERATIONS 

We now discuss some systems of experimental rela- 
vance to the physics studied above. Light alkali atoms 
that have been condensed are ''Li and ^^Na, although the 
latter is not really light when taking ratios with other in- 
teresting systems. However, recent advances in studies 
of metastable Helium, '^He* [s^, and the potential of 
those experiments to make mixtures of ^He* and ®^Rb 
make this a very interesting system [s^. Likewise, the 
mixture of '^Li and ^^^Cs [i^ or '^Li and ®''Rb [ilj would 
be relevant, and also potential rnixtixres of Ytterbium 
isotopes (mass numbers 168-176) |42', 43] and ^Li or ^He 
would be very favourable. 

The results we have presented here show that the Efi- 
mov spectrum can be modified at large distances by the 
presence of a condensate background, and also that there 
is no effect at short distance from the condensate. If we 
insert the definition of ^ into Eq. (|16p , we find the follow- 
ing formula for the expected number of three-body Efi- 
mov states at unitarity (which effectively assumes that 
|a| ^ ^) for our setup 



So 
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, (17) 



where hq is the condensate density, as is the scattering 
length of the condensed bosons, and rvdw is the two-body 
van der Waals length associated with the interatomic po- 
tential of the two heavy impurities. The unit ao is the 
Bohr radius. If we consider the ^He-^^Rb or ^Li-^'^'^Cs 
cases, the mass ratios are roughly the same and we have 
So = 1.98. For the ^Re*-^'^Rh case, ub ~ 142ao 113, HI 
and rvdw(*''Rb) ^ 83ao [1^, and we obtain Nb ~ 3.2 
at no — lO^'^cm"'^. One order of magnitude increase in 
no or aB brings this down to Nb ^ 2.5, so we see a 
sizable effect. Of course we are assuming that there are 
non-overlapping resonances in *He-^^Rb and ^^Rb-^^Rb 
which is currently unknown. The example of ^Li-^'^'^Cs 
is slightly more complicated since ''Li has attractive in- 
teractions at zero magnetic fields [s!]. However, let us 
for the moment assume that we can tune the scattering 
length away from the attractive region and also find a 
good resonance in the Li-Cs system (as recently done for 
the closely related case ^Li--'^'^'^Cs [44l |45|'). If we assume 
that qb ^ lOOflQ and use r^dw ~ lOlao [3a|, we find more 



or less exactly the same scenario as in the ^He-87Rb case. 
We thus see that the effect of condensation of light par- 
ticles in heavy-heavy-light three-body systems should be 
accessible in current experiments. 

Other mixtures are being pursued at the moment that 
are relevant for our purposes since they contain Bose 
components that can be condensed, but for which the 
mass ratios are too large for the Born-Oppenheimer ap- 
proximation to be fully justified. Examples are ^'^Na 
mixed with ^"K [Ifij and *'^Rb with ^^^Cs mixtures [47j . 
For the latter system the experiments that are probably 
most interesting for our purposes are of the few heavy 
impurities in a condensate of the lighter species by the 
Widera group [H, . While the approximations used 
here are not expected to be accurate in relation to these 
less mass imbalanced mixtures, we do expect that simi- 
lar signatures should occur, and that the Efimov effect, 
if present, will be modified by a condensate background 
when ^ is sufficiently small. 



VI. DISCUSSION 

We have considered the influence of a condensate back- 
ground on three-body bound state physics in the case of 
two heavy impurity atoms embedded in a condensate of 
light particles and assuming that the light-heavy inter- 
action is short-ranged. Using the Born-Oppenheimer ap- 
proximation we calculate the modification of the heavy- 
heavy interatomic potential when the light particle dy- 
namics is integrated out. Our results demonstrate that 
this potential is strongly modified at length scale corre- 
sponding to the condensate coherence length, and even- 
tually turns from attractive inverse square (necessary for 
the Efimov effect), through zero, and then most likely 
into a repulsive potential at very large distance. These 
findings indicate that the coherence length must be con- 
sidered as a length scale when estimating the potential 
for such systems to form universal three-body bound 
states. In the case where the heavy-light interaction is 
resonant (infinite scattering length), the coherence length 
replaces the scattering length in the famous Efimov for- 
mula, Eq. (ITBl) . We have estimated the effects of our 
findings on experimental mixtures of current interest and 
find that by tuning interaction strengths and condensate 
density, it should be possible to manipulate the number 
of universal three-body states. 

In future studies it is necessary to go beyond the Born- 
Oppenheimer approximation in order to not rely on very 
mass-imbalanced systems. Furthermore, two or three of 
the constituents may be condensed and it would interest- 
ing to study the effects that this will have on the three- 
body spectrum. From the current study we would expect 
that some combination of coherence lengths and scatter- 
ing lengths would decide the number of bound states. It 
would also be interesting to study the addition of degen- 
erate Fermi components. Here the coherence length is 
replaced by the Fermi momentum 0| ■ 
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Another important way in which to improve the cur- 
rent study is in the description of the Bose condensate 
through the Bogohubov formaUsm. If we have a strongly 
interacting condensate one needs to consider modifica- 
tions here. Also the neglect of states with quasi-particles 
(phonens) may not be justified then and we need to con- 
sider instead a starting point in line with the Frohlich 
polaron system JSQ]. However, we do expect to see qual- 
itatively the same physics, i.e. that presence of back- 
grounds will modify the universal three-body spectrum 
and in many cases suppress bound state formation. 

The essential assumption used here was that the light 



particle had a linear dispersion at low energy. This is 
reminiscent of the low-energy dispersion of a system like 
Graphene [5l| or the surface of a topological insulator [sll 
where it is electrons that have linear dispersion at low- 
energy. Adding impurities to such systems and studying 
whether bound states induced by the electronic surround- 
ings with linear dispersion is an interesting prospect. 

This work was supported by the Danish Agency for 
Science, Technology, and Innovation under the Danish 
Council for Independent Research - Natural Sciences. 
We are grateful to Manuel Valiente for feedback on the 
manuscript. 
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